Wigner crystallization in Na3Cu204 and NagCusOio chain compounds 
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We report the synthesis of novel edge-sharing chain systems Na3Cu204 and NagCusOio, which 
form insulating states with commensurate charge order. We identify these systems as one- 
dimensional Wigner lattices, where the charge order is determined by long-range Coulomb interac- 
tion and the number of holes in the d-shell of Cu. Our interpretation is supported by X-ray structure 
data as well as by an analysis of magnetic susceptibility and specific heat data. Remarkably, due to 
large second neighbor Cu-Cu hopping, these systems allow for a distinction between the (classical) 
Wigner lattice and the Akp charge-density wave of quantum mechanical origin. 

PACS numbers: 71.10.-w,71.28.+d, 75.10.Pq 
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The role of strong electron correlations and the con- 
comitant appearance of spatially modulated charge struc- 
tures constitutes a central issue in current solid state 
physics Q. The most prominent example are charge 
stripes in high-Tc superconductors, which have been dis- 
covered first as static modulations in the Cu02 planes 
of Lai.6-a:Ndo.4Sr2;CuO40. The high-Tc enigma has 
at the same time stimulated the fiorishing field of low- 
dimensional cuprates composed of Cu-0 chains and lad- 
ders Q with the hope to gain new insights overlooked so 
far in the layered cuprates. One important aspect is the 
role of long-range Coulomb interaction in strongly corre- 
lated systems, which gives rise to metal/insulator stripe 
structures in organic charge-transfer compounds^]. 

Following a recently discovered new route in the syn- 
thesis of alkalioxometallates [3, we were able to syn- 
thesize several members of a new class of quasi- ID 
cuprates Nai+a;Cu02. These intrinsically doped edge- 
sharing chain systems provide a unique opportunity to 
study the condensation of charge order (CO) at high 
temperature and the formation of spatially modulated 
Heisenberg spin systems at low temperature. Edge- 
sharing chains are also building blocks of the intensively 
studied Sri4_a;Caa;Cu2404i system; due to exchange of 
electrons with ladders, however, the degree of doping 
is difficult to determine j^. Here, we argue that these 
doped chains can be understood as realizations of one- 
dimensional Wigner lattices (WL) 0, as introduced by 
Hubbard ^ in the late 70 's in connection with TCNQ 
charge transfer salts. He suggested that the distribu- 
tion of electrons is controlled by the Coulomb interaction 
rather than by the kinetic energy ('^ band width), such 
that they form a generalized Wigner lattice on the under- 
lying TCNQ chain structure. This view suggests a strik- 
ingly different nature of charge excitations, namely as 
domain walls with fractional charge rather than particle- 
hole excitations as in common metals and semiconduc- 
tors. Hubbard's proposal, however, can be challenged on 
the grounds that the resulting periodicity of charge mod- 
ulation can alternatively be explained by a ikp charge 

density wave (CDW) I1ESE3 

arising from short-range 



interactions alone and an instability of the Fermi surface, 
where kp denotes the Fermi momentum. 

The edge-sharing arrangement of CUO4 squares meets 
the WL criterion of small band width in an optimal way 
due to the almost 90° Cu-O-Cu bonds (Fig.l). Unex- 
pected complexity is added because, apart from a small 
nearest-neighbor hopping matrix element ii, there is also 
a second neighbor hopping t2, which is larger as a con- 
sequence of the structure. While this unusual feature 
does not affect the classical WL order imposed by the 
Coulomb interaction, it changes the Fermi surface topol- 
ogy, and thereby allows to distinguish the WL from the 
CDW on the basis of the modulation period. 

These systems provide a first example where an un- 
ambiguous distinction between the generalized WL and 
a Fermi surface related ikp CDW is possible. We also 
show that for these edge-sharing compounds even the 
magnetic and thermodynamic properties can only be ex- 
plained by invoking a WL ground state emerging from 
the truly long-ranged Coulomb interaction. 

Samples were prepared by the azide/nitrate route j^. 
As a source for the alkalimetal component, mixtures of 
the respective alkali azides and nitrates (or nitrites) are 
used instead of the alkalioxides. Conveniently, besides 
the metals' ratio, also the oxygen content and thus the 
degree and kind of doping of the desired product can 
be effectively controlled using the weighed portions of 
the starting materials. Following this procedure the ti- 
tle compounds Na3Cu204 {x = 1/2) and NagCusOio 
{x = 3/5) have been prepared as microcrystalline, pure 
phases in gram-amounts. The new oxocuprates (H/HI) 
belong to the compositional series Nai+a;Cu02, with the 
end members NaCu02 ^3 ^^'^ the still elusive Na2Cu02. 
The most prominent structural feature, common to all 
representatives known thus far, is a one-dimensional 
polyanion CUO2 ~ constituted of Cu04-square units shar- 
ing edges in trans-position. These anionic entities are em- 
bedded by sodium ions which achieve coordination num- 
bers of 4-6 with Na-O bond lengths ranging from 227 
to 279 pm. The geometric data as determined by sin- 
gle crystal structural analyses give clear evidence for a 
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FIG. 1: Structure of edge-sharing copper-oxygen chains along 
6-direction in Na3Cu204 (a) and NagCusOio (b), where Cu'^^ 
and Cu^^ are marked by green and red circles, and (Na^) 
by large blue (grey) circles, respectively, (c) View on the a — c 
plane of Na3Cu204. For convenience we introduce a cartesian 
system a' ,b'\\b,c' in addition to crystallographic coordinates. 



CO at the Cu sites (Fig. 1). Based on the Cu to O 
bond lengths, one can unambiguously distinguish Cu'^+ 
and Cu^+ sites. The way of linking the primary struc- 
tural units together with the variations of the copper to 
oxygen distances inevitably leads to deviations of the O- 
Cu-0 angles from the ideal 90°. As monitored by differ- 
ential scanning calorimetry measurements CO disappears 
above the WL melting temperature Tm = 455 and 540 
K for Na3Cu204 and NasCusOio, respectively. DC and 
AC conductivity measurements show also a clear transi- 
tion from an Arrhenius behaviour below Tm to an almost 
temperature independent conductivity regime above T,„. 

For a theoretical analysis one has to recognize that 
Cu^"*" is in a configuration with spin 1/2, while Cu'^^ 
is in a d^-ligand hole {(PLh) state, also known as Zhang- 
Rice singlet state |0. In contrast to high-Tc cuprates 
the edge-sharing geometry (Fig.l(a,b)) leads to strongly 
reduced hopping matrix elements. This sets the stage for 
the long-range Coulomb force as dominant interaction 
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where the on-site interaction U suppresses charge fluctu- 
ations involving Cu^+ {d^'^) configurations. Here we asso- 
ciate the (PLh (d^,d^°) ionization state with (1,2) elec- 
trons, respectively, and rii^cr (f =Tji) counts the num- 
ber of electrons with spin a, while rii — n^^i + 
The Coulomb interaction Vi in general is screened by 
the polarization of neighbouring chains as well as by 
core electrons Jsj. Here we assume a generic Coulomb 
law Vi = ^ — 1,2,....^^. Crucial for the following 
is that the interaction is long ranged and convex, i.e., 
V/' = Vi-i-2Vi + Vi+i>0. 
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FIG. 2: (color online) Wigner charge order resulting from 
Coulomb repulsion and associated modulated Heisenberg spin 
structure for x — 1/2,3/5,5/8 and 3/4 doping (a-d). The 
spin-1/2 of Cu^^ (arrows) is responsible for magnetism; Cu''"'" 
(circles) is nonmagnetic. The spin arrangement is that ex- 
pected for ferromagnetic exchange Ji and antiferromagnetic 
J2, where the two excitations dP -(f d?Lh ■ d^*^ contributing 
to J2 are indicated by dashed arrows in (b). The charge unit 
cells (shaded) contain 2, 5, 8, 4 sites, respectively. The struc- 
tures (a) and (b) are realized in Na3Cu204 and NagCusOio, 
respectively, and the dashed circles in (a*) and (b*) indicate 
charge excitations (x V in these structures, (d) Typical mod- 
ulation for charge stripes in cuprates at doping 1/8. 



For commensurate doping concentration x = rnjn 
the interaction Vi selects a particular CO patternQ. 
This pattern is immediately obvious for filling fractions 
X = 1/2 and 3/4 (Fig. 2(a,d)) which involve an equidis- 
tant arrangement of the Cu'^"'" sites (red circles in Fig. 
2). For a general ratio x = m/n this leads to complex 
structures with unit cell size n (in units of the Cu-Cu 
distance b' = 1). In case of a; = 3/5 we encounter in Fig. 
2(b) the charge order observed for NagCusOio. Charge 
localization, however, is not perfect in Wigner insulators 
as electrons still undergo virtual transitions to neighbor- 
ing sites (Fig. 2(a*,b*)) in order to retain partially their 
kinetic energy. The energy of the lowest excitations and 
the impact of kinetic energy depend strongly on x, e.g., 
the energy of the excitation in Fig. 2(a*) relative to the 
ground state Fig. 2(a) is ^ V2 while the excitation for 
a; = 3/5 in Fig. 2(b*) is ~ V^', about an order of magni- 
tude smaller. To investigate the role of kinetic energy we 
explore the dynamics of electrons starting from the ID 
Hubbard- Wigner model Hhw = Hcoui + Hxin H , where 
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describes the hopping of an electron with spin a. Due to 
the almost 90° Cu-O-Cu angle the hopping ti between 
nearest neighbor Cu sites results mainly from direct d — d 
exchange, while t2 originates from hopping via a Cu-0- 
0-Cu path 15] (Fig. 1), leading to the remarkable fact 
1^2! > We adopt here as typical values ti ~ 63 
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FIG. 3: (color online) (a) Static charge structure factor Sc{q) 
calculated for interacting spinless fermions and ti — 0.04^, 
t2 = 0.061^ at x = 0.6 indicates the instability towards WL 
ordering with modulation qw~l-2n (solid line as guide to the 
eye). Inset shows the temperature dependence of Sc{<lw)- (b) 
Comparison of electron dispersions for a system with nearest 
neighbor hopping ti{— 1) only (dashed line) and for ti = 
1, t2 — 1.5 (solid line) with the Fermi energy indicated by 
a horizontal line, (c) While the charge susceptibility Xcil) 
for noninteracting spinless fermions in the ti model shows a 
logarithmic divergence at 2ksF ~ q-w, the singularities of the 
ti-t2 model are at different momenta. 



meV, t2 ^ 94 meV, derived from ab-initio band structure 
calculations for the Cu^+ edge-sharing reference system 
Li2Cu02 These values are indeed much smaller than 
our estimates for U ~ 3.8 eV and V ^ 1.5 eV based on 
optical data for Li2Cu02 [isj . 

Using exact diagonalization we have calculated the 
static structure factor Sc{q) for chains up to 25 sites. 
The peaks of Sc{q) at — \.2ti (and at 0.87r=27r-g^) 
are characteristic for a WL modulation at x = 3/5 (Fig. 
3(a)). Since the magnetic energy scale is much smaller 
than the Coulomb interaction V ^ the charge structure 
can be determined by disregarding spin degrees of free- 
dom, namely in terms of spinless fermions (SF). The SF- 
CDW arises as an instability due to low energy scat- 
tering between the two Fermi points at iksF which 
lead to a singularity in the charge susceptibility Xc (q) 
q = 2ksp{=AkF)- Subsequent inclusion of interactions 
is expected to change the character of a singularity, yet 
not the momentum at which the singularity occurs. This 
explains the origin of the Akp CDW instability in inter- 
acting ID systems and, as Akp — 27rm/n, the coincidence 
with the modulation period n of the WL. 

The equivalence of WL and CDW periodicity disap- 
pears when t2>at\ (a — 0.38 at x = 0.6); then there are 
four instead of two Fermi points leading to new singular- 
ities in Xc(9) (Fig.3(b,c)) and to a shift of the original 
singularity away from Akp = qw In this case the stan- 
dard tools of many body theory such as bosonization [ll| 
suggest a change of the modulation period. This, how- 



ever, is not reflected in the experiment and is also not ob- 
served in our numerical results for Sdq) (Fig. 3(a)). In 
fact, the singularity remains at qw up to a value t2 ~ 4ti, 
a striking manifestation that the structure is robust and 
controlled by the long range Coulomb interaction [l7j . 
The calculated Sc{q) shows besides the main peak at q^ 
a weak higher harmonic feature aX q — 2g^ — 27r, con- 
sistent with the relative intensities of structural reflexes 
calculated from the experimentally determined Cu ion 
positions (vertical bars in Fig. 3(a)). Further results (in- 
set Fig. 3(a)) reveal that WL correlations persist up to 
temperatures fcsT 0.05^, consistent with experiment. 

The magnetic susceptibility x(T) was measured in the 
temperature range from 5 K to 350 K using a SQUID 
magnetometer. The data (Fig. 4(a)) reveal strikingly dif- 
ferent temperature dependences for the two compounds: 
x(r) for Na3Cu204 displays some similarity with a near- 
est neighbor Heisenberg antiferromagnetic chain [l8| . 
whereas the NagCusOio data show continuous increase 
down to low temperature until its maximum near 25 K is 
reached. Both systems reveal antiferromagnetic correla- 
tions, yet magnetic order is observed only for NagCusOio 
at Tjv = 23.5 K. For a calculation of x{T) we assume that 
the spins remain fixed at their positions Ri as given by 
the structural analysis and by Hhw (Fig. 2). This leads 
to a generalized Heisenberg model 

Hhc^s ^\Y1 ^(^» - ^j)^' • (3) 

where the exchange constants depend on the distance 
|Ri — Rj l between the spins and on the direction (paral- 
lel to the chains or perpendicular). For x = 0.5 only the 
exchange constants J2, Ji, ■■■■ along the chains contribute, 
while for x = 0.6 Ji, J2, Jz, ■■■■ are relevant. Apart from 
the modulated spin pattern, superexchange [13 in WL's 
shows further novel features, namely fluctuations (a) of 
spin positions and (b) of antiferromagnetic exchange in- 
tegrals Ji ~ 'itf/{U +A{V)) due to the low energy charge 
fluctuations. The energy shifts A{V) due to Vi depend 
on the WL structure, and are in general different for 
left (right) scattering processes (see Fig. 2(b)). The ex- 
change constants may be estimated from the parameters 
specified above. The largest couphng is J2 ~ 100 — 200 
K resulting from superexchange via a Cu-O-O-Cu path. 
Because of the almost 90° Cu-O-Cu bond angle (13) Ji 
exchange is smaller. Yet in the case of Ji there are ad- 
ditional exchange contributions, which have to be taken 
into account. The most important of these involve O 

configurations with Hund interaction on oxygen |2nl |. 
such that the total Ji may become negative. 

We have used finite temperature diagonalization 
(FTD) ^ to calculate xiT) for Na3Cu204 studying 
chains up to = 48 sites (24 spins). Good agreement 
with experiment was achieved with J2 — 172, J4 = 17 K 
and 1.97 for the g- factor. The small antiferromagnetic in- 
terchain coupling J' = 34 K (Fig. 4(b)) is frustrating as 
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FIG. 4: (color online) (a) Temperature dependence of suscep- 
tibility xiT) for Na3Cu204 and x/3 for NagCusOio. Mul- 
tiplication by 1/3 in the latter case provides a comparison 
per spin. Theoretical results calculated by FTD are indicated 
by lines (exchange constants are given in the text). Sketch 
of interchain interactions used in the calculation of x{T)'- 
For Na3Cu204 (b) the magnetic exchange J' is frustrated 
due to triangular coordination, while it is not frustrated for 
NagCusOio both along a'-direction (c) and c'-direction (d). 



at high temperature x(T) is controlled by the individ- 
ual spins, at low temperatures only effective spins 1/2, 
formed by three spins which are coupled by the largest 
interaction J2, are relevant (Fig. 2(b)). The compos- 
ite nature of the effective spins and their internal energy 
level structure becomes relevant in the intermediate tem- 
perature region. The interaction between these effective 
spins is renormalized to J^'^'^ = ^[Ji — J3 + jJd]- Thus, 
the peculiar low-T behavior of x(T) is controlled by an 
interplay of the small Ji'^^ between effective spins, re- 
sulting from WL order, and the interchain interactions. 

Finally, the new compounds provide a unique oppor- 
tunity to study the competition between two entirely dif- 
ferent states, the classical WL dictated by the long-range 
Coulomb interaction and the CDW of quantum mechan- 
ical origin, i.e., resulting from a Fermi surface instability. 
These materials highlight the importance of long-range 
Coulomb interaction in strongly correlated systems, — 
and provide a one-dimensional test ground for the study 
of charge stripe formation. 

We thank E. Briicher and G. Siegle for the susceptibil- 
ity and specific heat measurements and L. Capogna, B. 
Keimer, W. Metzner and A.M. Oles for discussions. 



the correlations on b' chains are antiferromagnetic. The 
interaction J' has been taken into account in a mean- 
field approach^ xiT) = Xs{T)/il + z'\J'\xsiT)), where 
XsiT) is the result for a single chain, while z' counts the 
number of spins on neighboring chains (z' = 4). The 
frustrated geometry implies a spin gap Ag For the 
present parameter set we obtain ~ 4 K, consistent 
with our analysis of the specific heat data. A thorough 
analysis of the Curie contribution, using specific heat 
and susceptibility data independently, reveals that there 
are about ^1% impurities per Cu ion present in both 
Na3Cu204 and NagCugOio 23]. 

In the case of NagCusOio the correlations along the 6- 
chains turn out to be of ferrimagnetic nature. A full FTD 
calculation of a double chain was necessary to account 
for the AF interchain exchange J' along c'-direction (Fig. 
4(d)), while the coupling of chains along a'-direction (Fig. 
4(c)) was treated within mean-field. This implies fer- 
rimagnetic a'-b' planes, which are antiferromagnetically 
coupled along c'-direction. A satisfactory description was 
found for Ji = -22 K, J2 = 113 K, J' = 21 K and 
g = 2.0. The ferromagnetic interaction Ji is the origin of 
the strong enhancement of x(T) in the x — 0.6 system. 
The reduced value of J2 can be explained by low-energy 
charge fluctuations (Fig. 2(b*)), which involve a fluctu- 
ation of spin positions, i.e., J2^^ ~ (1 — p)J2 + pJi and 
a larger probability p than for x = 1/2. 

The strong variation of x{T) in NagCusOio should be 
interpreted as a crossover between two regimes. While 
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